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Abstract 

. This paper is concerned with the linear ODE in the form y'{t) ~ Xp{t)y{t) + b{t), A < and which represents 

' a simplified model of storage of the carbon in the soil. In the first part, we show that, for a periodic function 

ly^ . p{t), a linear drift in the coefficient b{t) involves a linear drift for the solution of this ODE. In the second part, 

CO ' we give sufficient conditions on the coefficients to ensure the existence of an unique periodic solution of this 

. differential equation. Numerical examples are given. 

o ; 

oo . 

^ Keywords: Ordinary differential equations, T-periodic function, linear drift, Cauchy sequence, series' estimates, 

• • . AMS subject classification: 34E05, 40A05. 



1 Introduction 



A lot of phenomena of evolution are described using ordinary differential equations ODE or systems in which 
the coefficients and/or the source terms are periodic. Let us mention some applications in physics (e.g. the 
harmonic oscillator , the resonance phenomena due to oscillatory source terms) , electricity (let us mention the 
famous RLC circuit with an oscillatory generator), in biology (circadian cycle), in agricultural studies (due to 
seasonal effects). 

The main question which is addressed here is when or under which conditions a slow perturbation of the 
coefficient in the ODE will induce a similar behavior on the solutions, in large time. More precisely, we are 
looking for conditions to ensure that a linear drift in the (periodic) coefficients of the ODE will lead to a 
linear deviation (and thus unbounded) of the solutions or, on the contrary, what kind of perturbation in the 
coefficients are compatible with a stable (bounded, periodic in large time) solutions. 

Although these questions can be applied to large number of applications, our original motivation concerned 
with the effect of climate change on the seasonal variations in the organic carbon contained in the soil as claimed 
at the end of the conclusion of Martin et al. [13]. 
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Since the readers may be not familiar with this domain. Let us recall some basis about the issues of the 
soil organic carbon (SOC) modelling. 

The spreading is one of opportunities for organic materials of human origin (sludge of filtering treatment 
station and derivative products) and agricultural (manure). The organic matter spread contain significant 
amounts of organic carbon which, after application, a fraction is permanently stored in the soil. ; The remainder 
is returned to the atmosphere as C02. The spreading can be accessed through the storage of organic carbon 
in the soil, helping to reduce C02 emissions (a major greenhouse gas effect) compared for example to the 
incineration of organic matter that returns all carbon into the atmosphere. The optimization of spreading of 
organic materials is important in reducing emissions of greenhouse gas effect. 

The dynamics of carbon in the soil, which determines the amount of organic carbon stored in soil and 
returned in the form of C02, depends on soil type, agricultural practices, climate and quantities spread. The 
soil organic carbon (SOC) plays an important role in several environmental and land management issues. One 
of the most important issues is the role that SOC plays as part of the terrestrial carbon and might play as 
a regulation of the atmospheric C02. Many factors are likely, in a near future, to modify the SOC content, 
including changes in agricultural practices [3,18,2] and global climate changes [9,5,8,10,11]. 

Understanding SOC as a function of soil characteristics, agricultural management and climatic conditions 
is therefore crucial and many models have been developed in this perspective. These models arc used in a 
variety of ways and after for long term studies [6,15,16,17]. The behavior of the SOC system, over a long term 
and assuming that the environment of the system (inputs of organic carbon, climatic conditions) is stable, is 
reported to tend towards steady or periodic state. 

Some of he SOC models have been formulated mathematically [14,4,1,12,13]. Here we consider the RothC 
model [7,13] which consists in splitting the soil carbon into four active compartments. Under a continuous 
form it can be written as 

^=p{t)AC{t)+B{t), (1.1) 

where C (t) is a vector with 4 components, each corresponding to a compartment storage of carbon in the soil: 
DPM (decomposable plant material), RPM ( resistant plant material), BIO (microbial community) and HUM 
( humus). These components indicate the amount of carbon stored at the moment t. 
In (1.1) 

i) p{t) is a function indicating the speed of mineralization of soil, which results in C02 emissions. 

ii) A is a matrix that can represent the percentage of clay in the soil and speeds of mineralization in each 
compartment of C in the soil. 

iii) B(t) indicates the amount of carbon brought in the soil (amount spread per unit time). 

Let us recall that the initial goal of this study was to understand how long term evolution on climatic 
data imply variation in large time on the solution of ODE with periodic coefficient and/or source terms with 
a (linear) drift. The mathematical tools involved in this paper are rather classical and simple but there are, 
up to our knowledge, very few literature on the subject. 

In the reminder of this paper, we shall consider a simplified case with a single (scalar) equation but the 
extension to a diagonahzable system (as it is the case for the A matrix in [13], see eq (5)). 

The main results in this paper are, first (section 2), to prove that a linear drift in the coefficients leads to 
a linear drift in the solution over a large time behavior. Second (section 3), we try to find sufficient conditions 
on the coefficients in order to get the existence of an unique periodic solution. These results are illustrated by 
numerical tests in section 4. 

The extension to linear differential system is straightforward (by diagonalizing the system). The extension 
to partial differential equation arc under study. For example, in the case of the linear heat equation, using 
Fourier Series, one can expect the same kind of results. 
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2 Asymptotic behavior of the solution 

In this paper, we consider the linear differential equation 

?/ (t) = \p [i] yit) + b{t), 0<t< +00, (2.1) 

where A < 0, p{t) and b{t) are given functions satisfying the following conditions 
(Ai) p(t) is an T-periodic function, with T > fixed. 

{A2) there exist the T-periodic function, (3{t), such that 

b{t + T) = b(t)+p{t), Vte[0,oo). (2.2) 
The general solution of (2.1) has the form 

y{t) = e«W jci + J* 6(s)e-»(^)ds| , (2.3) 

where Ci is a constant and 

a{t) A / p{s)ds. (2.4) 







In this section, we prove that there exists a unique solution yoo{t) of (2.1) satisfying 

yoo{t + T) = y^{t)+j{t),yt€ [0,00), (2.5) 

where 7(4) is an T-periodic function. Let first remark that 

Lemma 2.1. Let b : i?+ — > R. The following properties are equivalent : 
(a) 3l3{t) periodic with period T such that b(t + T) = h{t) + /3(t), Vt > 
(h) 3b{t) periodic with period T such that b{t) = b{t) + j^f3{t),yt > 

Proof. Let us first prove (a)=^'(b). Choose b{t) = b{t) — ^f3{t). Then b{t) is periodic with period T since 
b{t + T)= b{t + T) - i^^{t + T) = b{t) + f3{t) - (^ + l)/3(i) = h{t) - ^/3(<) = b{t). 

Conversely (b)^(a). Using h{t) = b{t) + ^f3{t), we have b{t + T) = b{t + T) + ^^^p{t + T) = 
b{t) + (7^ + l)P{t) = b{t) + ^p{t) + I3{t) = b{t) + (3{t). This concludes the proof. □ 

Let us now prove the main result of this section. First, we state the some lemmas 

Lemma 2.2. Let (Ai) hold. Then 

a{t + nT) = a{t) + a{nT) = a{t) + na{T),\lt > 0, n e N. (2.6) 

Proof. From (2.4) we deduce that 

ait + nT) = X / p{s)ds + A / p{s)ds = a{nT) + A / p{s)ds. (2.7) 

Jo JnT JnT 

On the other hand, by the assumption (Ai), we have 

pt+nT pt pt 

A / p{s)ds = A / p{s + nT)ds = X p{s)ds = a{t), (2.8) 

JnT Jo Jo 

and 

n-l p(k+l)T n-1 pT 

a{nT) = A ^ / p{s)ds = A ^ / p{s)ds = na{T). 



(2.9) 



fe=o " «^ fe=o " 

Combining (2.7)-(2.9) we have (2.6). □ 



3 



Lemma 2.3. Let assumptions (Ai), {A2) hold. For n G Z+ and t G [0, 00), we put 



(pt+nT ^ 
2/(0) + 1 b{s)e-'^^^Us 



Then, 

yn{t)=yoc{t) + 6n{T) + n^{t), 

where 

{po(T) fT a(T) pT r-t 

{rMT) f-T a(T) rT ) 

y{0) + —-^ 5(s)e-«Wrfs+- /3(5)e-«Wrf4, 
i — e ^ ^ Jo (1 - e'^v-' )j Jo J 

{pa(T) pT i-t ] 

^— ^ (3{s)e-^(^Us + l3{s)e-<^Us \ . 

Proof. By the assumption [A^), it follows from (2.10) and the lemma 2.2 that 

y^{t) = e"(*) |y„(0) + ^* 6(s)e-'^(^)ds + n^* /3(s)e-'^(^)ds 
On the other hand, we have 



and 



2/„ (0) = y {nT) = e'^("^) <| y (0) + / h [s) e-'^^'Us 







= e""(^) i 2/ (0) + V / & (s + yfcT) e-'^(^+'=^)ds 
I fc=o-^o 

'(0)+ / 6(s)e-"('^)cis^e-'=«(^) + / /3 (s) e-"(")ds ^ fee 

Jo u_a Jo u-a 



fc=0 • 

T n — 1 n—1 



fe=o fe=o 

By using the following equalities 

1 „-na(T) 

g-fea(T) _ 

1 - e-<T) ' 

fe=0 



and 



e-a(T) e-(n+l)a(T) g-na(T) 



fc=o (e-«(T)_i)2 (e-a(T)_i)2 e-«m-l' 

thus we obtain 

P<T) fT a(T) rT 

f fT a{T) fT 

+ e-(-)|.(0)-^-^y^ 6(.)e-«(«).a + ^^--^y^ Hs)e-^^^<is ^ 



e 

■ n 



a(T) i-T 

(3{s) e-^'^ds 



'l-e"(^) 7o 

Combining previous equalities, we obtain (2.11). The proof of Lemma is complete. 



Now, we state the main theorem 
Theorem 2.4. Let (Ai), {A2) hold. Then, there exists a unique solution yoc{t) of (2.1) such that 

yoo{t + T)=y^{t)+-f{t),yt€ [0,00), 
where 7(f) is an T-periodic function, defined by 



j{t) = e'^(*) 
Proof. For n G Z+ and t > 0, let us define 



pa(T) rT ~) 

Y:r^l /3(s)e-»(«)ds + p{s)e-<^^ds\ 



def 

Un{t) = yn{t) -n-f{t). 



It follows from (2.11)-(2.14) and (2.17) that 

lim Un {t) = yoo {t) , yt e [0, +00) . 

n — ^+00 

It is clear that yoc{t) is a solution of equation (2.1) satisfies the initial value 



By (2.18), we have 



yoo {t + T)= Um Un {t + T)= Hm {y„ {t + T)- 717 {t + T)} . 

= lim {y„+i(t)-(n + l)7(t)+n[7(t)-7(t + T)]+7(t)} 

n — ^+00 



On the other hand, by the periodicity of f3{t), we get 
7(i + r)=e«(*+^) 



ea(t) 



a(T) rT nt+T \ 

-—^^j^l3{s)e-<^Us + j^ p{s)e-'^(^Us\ 
^— ^ /J(s)e-«Wds + /3(s)e-«Wrf5U7(t). 



Combining (2.20) and (2.21), we obtain 

2/oo(i + T)= lim M„+i(t)+7(t) = 2/oo(t)+7(i)- 

n— >+oo 

Uniqueness 

Now, let ^(f) be the solution of (2.1) corresponding to the initial value y (0) = A and 

y(t + T) = y(t)+7(t), 

where 7 (i) is an T-periodic function. Then y* (t) = y^o (t) — y {t) satisfy 

y' (t) =Xp(t)y(t), 0<t< +00, 
y{0)=L (T) - A, 



and 

y* it + T) = y* (t) + 7* (t) , 7* + T) = 7* (t) , Vt > 0. (2.25) 

It follows from (2.24) that 

y* (t) = {L (T) - A) e"(*) ,\/t>0. (2.26) 
From (2.25) and (2.26) we deduce that 

7* {t) = -{L (T) - ^) (1 - e''^^)) e"W, > 0. (2.27) 

Combining (2.25), (2.27) we get A — L{T). By the uniqueness of Cauchy problem, the proof of theorem 2.8 is 
complete. □ 

Remark If we consider the equation (2.1) where the assumptions {Ai) and (^2) are replaced by 
{A\) b {t) is an T-periodic function, with T > fixed 

{A'2) there exist a T-periodic function, a{t), such that 

p{t + T) = p{t)+a{t), VtG[0,oo). (2.28) 

In that case it is clear that there does not exist a solution which has the same property as the function p{t), 
for instance if we consider the example with b{t) = 0. Here the solution of (2.1) tends to as f — > +00. 

3 Sufficient conditions for the existence of periodic solutions 

In this section, we consider the following assumptions for the functions p{t) and b{t) of the ODE (2.1) 

(Bi) p € L;^^^(0,oo), p{t) > and there exists a constant T > such that 

p{t + iT) = p{t) + ai{t), e [0,+oo),Vn e N, (3.1) 

{B2) Assume that the function t ^ h(t)e-^^« belongs to Ll^^{i), 00) and 

b{t + iT) = b{t)+ I3i{t), Vie [0,+oo),VneN, (3.2) 

where the functions sequences {Q!i(i)}, {/3i(i)} satisfy some conditions specified later. Note that, by (3.1), (3.2) 
we also deduce that G L]^^{Q,oo) and /3i(-)e--^/o g L]^^{Q,oo). 

Let us define the function sequence {y„(i)}^o by 

yn{t)=y{t + nT), (3.3) 

where y{t) is a solution of equation (2.1) corresponding to the initial value y(0) = C\ and defined by (2.3). It 
can be proved that if p{t) and b{t) are the T-periodic functions then the function sequence {?/„(i)}^o converge 
to 2/00 (i) where 

{pa{T) fT ft ] 

Y^r^l 6(s)e-«(^)rfs + y^ 6(s)e-«(^)rfsi, (3.4) 

2/oo(i) being a unique T-periodic solution of equation (2.1). Now, we will extend this result to the case where 
p{t) and b{t) satisfy the conditions {Bi), (-B2), respectively. 
Let us now state some lemmas. 
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Lemma 3.1. Let (Bi) hold. Then we have 

a{t + nT) = a{nT) + a{t) + X [ a„(s)ds, Vn € Z+, Vi > 0, 

and 

n — l 



a(nT) = na(r) + A / ak{s)ds. 



Proof. From (3.1) and (2.4), it follows that 

pnT pt+nT 

a{t + nT) = A / p{s)ds + A / p{s)ds 

Jo JnT 

= a{nT) + X J p{s + nT)ds — a{nT) + a{t) + X I a„(s)ds. 
Hence, we obtain (3.5). On the other hand, from (3.5) we can deduce that 

a {{n + 1)T) ^ a{nT) + a{T) + A / a„(s)ds. 



By recursion, we obtain p.6p from 
Lemma 3.2. Let (Si), (B2) hold. Then 

2;„(t)=e''W+^""(*){y„(0) + ^„(0}, 

where ^ ^ 

a,(i)=/ a,{s)ds, and A(0 = / (bis) + P,{s)) e-''^'^'^"'^'Us, i&N. 



Proof. From (3.3), (2.3) and (3.5), we deduce that 

f nnT+t 



(•riT rnT+t j 

„a(nT)+a(t)+A/>„(.)ds I _^ / t{s)e-''^''Us + / &(s)e-"(")ds 



nT 



•nT 

Ci+ / &(s)e-''("'ds 



<.nT+t 
JnT 



I [ ^ ' -Jo 

( /.riT+t ^ 

= ga(t)+A/„ a„(.)ds ) y^^(o) ^ ga(nT) 6(5)6-^(^)^5 I . 

On the other hand, by change variable, it follows from the assumptions {B2) and (3.5) that 

nT+t (.t 

b{s)e-''^''>ds^ / 6(T + nr)e-"(^+"'^)dT 

nT ^0 
/•t 

(fe(T) +/3„(r)) e-''("^)-'^M-^-/'o' dr. 
Combining p. lip and p.l2p . we obtain (|3.9p - (|3.10p . The proof of lemma 3.6 is complete. 



Now, replace t = T in p.9p . we get 

yn+l(0) = yn{T) = fJ. n {yn{0)+dn), (3.13) 

where 

M„ = e"(^)+^"", (3.14) 

and 

eeS,(T), A=A(T), VigN. (3.15) 
By the recurrent, from (|3.13p . we obtain the following corollary 
Corollary 3.3. Let (Bi), {B2) hold. Then 

n 

2/„+i(0) = 5]5i"^^„_,+i +<5(")y(r), (3.16) 

where 

^(n) ^g,a(T)+A(3„+a„_i + ...+a„ + i_,)^y^^ ^^2,...,n. (3.17) 

Now, to obtain the convergence of sequence {yn{0)}^=oy we make the following assumptions 
Assume that {ai{t)}°^^ and {Pi{t)}^i satisfies 

oo 

< a^+i{t) < a^{t),\/t>0 and ^ e^^^^^^Sj < +00, (3.18) 

1=1 

{B3) 

00 

< P^+lit) < ^3^{t),yt>0 and ^e-"'^)^ < +00. (3.19) 

i=l 

Remark 3.4. By the assumption {B3) we deduce that 

lim an{t) = 0,Vi > 0, (3.20) 

n — *oo 

and ^ 

lim X{t) = I h{s)e-''^'^ds.\/t > 0. (3.21) 

Indeed, the proof of (|3.20[) is straightforward and we omit it. From (|3.10[) we have 

X(t) - ^ &(s)e-'^(^)ds = ^ (6(s) + Pnis)) e""'^'') |e"^""(^) - l} + ^ /3„(s)e"'^(^)ds. (3.22) 
Using the following inequality 



|e^ - l| < |X|el^l, VXeR, (3.23) 
it follows from ([XTg]) , (|XTg)) and ([X^ that 



^„(t) - / &(s)e-'^(^)ds 



<||A|e-^"^y [|&(s)|+/3i(s)]e-'^(^)dsja„(t)+y /3„(s)e-'^(^)ds. (3.24) 
Hence, we obtain ((3?2T|l from (|3?23| by ((3?T9)) and (|3?20l) . 

Lemma 3.5. Lef (-Bi) - (S3) /loW. Then, the sequence {2/n(0)},j5°^o ^■^ convergence and we have 

lim y„+i(0) = -^—^ r b{s)e-'^^^^ds ^ L{T). (3.25) 
rj— >+oo 1 — e"'--' Jq 



Proof. From the assumption {B3) and (|3.17p . we have 



when n —^ 00. Hence, 



n / p(s)ds + Si + 32 
Jo 

hm (5i"^2/(T) = 0. 



Now, put 



We shah prove that {Z„} is a Cauchy sequence. Let m, n G Z+, m > n, then 



m n 



'm— 2+1 



'i— n+1 



It is clear that 



^(fc) ^ gm(T)_gA(3,+a,_i+...3,_,+i) < VfceN, i = l,2,...A:, 
with pi = e"^'^) e (0, 1) and by (|XT5)) - ([XTg)) . we have 

< / {\b{s)\+f3i{s))e-''^'^-^^'^''>ds = K{T), G N. 



Estimate Ji. 

It follows from (^30)) and (P3T|) . that 



z— n+l 



m — 2+1 



< if (T) E P\ ^ 0' 



z— n+l 



when m, n ^ +cxo, by X^i^i Pi < +00. 
Estimate J2. 

By p.30p - p.3ip . we can estimate J2 as follows 



1=1 

n 

< E ^1 /5m-i+l 



/3m-i+l — Pn-i+l 
Pn-i+1 



2^1 

n 



(™) _ ^(") 



/3n- 



■i+1 



(m) ^(n) 



On the other hand, it follows from (|3.10p that 



g Aa,„_i+i(s) _ g-A3„_i+i(s; 



Using the following inequality 



\ < \X-Y\,yX,Y > 0, 



we deduce from (|3J8| . ((3l9)) and (lOI)) that 



< |A| / {\bis)\ + e-"(«) (a„_,+i(s) - ds 



(3.36) 



<|A|jy^ (|6(s)|+/?i(s))e-'^(^)ds 



Hence, 



(3.37) 



On the other hand, we have 



Ei^ n+l -i^ ^ n+1 -i^ A 



(3.38) 



and 



i=i -^0 ,,^1 Jo Jo 

By the assumption (i?3), it foUows from p.37p - p.39p that 

0, when m, >oo. (3.40) 



Pi |/3m-i+l — /3n-i+l 
i=l 



Now, we shall estimate the term 



-(rn) An) 



It follows from p.lTp and the inequality (|3.35p that 



(3.41) 



> j^n — 1+1 j—m—i-\-l 



j—n—i-\-l 



Therefore, we deduce from (|3.4ip that 



E " ^1"^ ^ 1^1 E p\ E = 1^1 E E pi 



n-j+l 



(3.42) 



2—1 j— n— 2+1 



|A| 



n+l 
Pi 



n+l oo_ 
Pi 



I -Pi 



^Si — Q-i I < |A| — P]^ *Si ^ 0, when m, n ^ oo. 



Combining (|3.29p . p.32p . p.33p . (I3.40p and p.42p . we deduce that {Zn] is Cauchy sequence. Hence, it is clear 
that there exists the lim„^oo yn+i(0), by p.l6p and p.26p . Passing to the limit in p.l3p by p.20p - p.2ip . we 
have ^ 

hm j;„+i(0) = Y^^^ / b{s)e-<^^ds = L{T). (3.43) 



n — >+C30 

The lemma 3.5 is completely proved. 



□ 
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Theorem 3.6. Let {Bi)-{B3) hold. Then there exists a unique T-periodic solution of equation (2.1), yoa{t), 
defined by 

yoo(t) = e<'^ \l{T) + / &(s)e-"(^)dsl . (3.44) 

"'0 

Proof. First, by passing to the limit in (1^ . it follows from ([X^ . ((H?^ and ([5?^ that 

lim yn{t) = yoo{t),yt > 0. 

n — ^oo 

It is clear that yoo{t) is a solution of equation (2.1) and satisfies the initial value ?/(0) = L{T). Moreover, since 
y^{t + T)= lim yn{t + T)= lim y{t + T + nT) = lim yn+i{t) = yoo(^), Vi £ [0, +oo), 

n — >oo n — *oo n — 'oo 

we deduce that yoo(i) is a T-periodic function. We can also show as in section 2 that there is uniqueness of 
such a solution yoo{t). 

□ 



4 Numerical results 

In this section the following Cauchy problem (2.1) with the following choice 

A -1, yo = 1, p{t) ^ sin^t, b{t) = t 



(6.1) 



Clearly, p{t) and b{t) satisfy p{t + n) = p{t) and b{t + tt) — b{t) + tt i.e. the assumptions of section 2 with 
T — TT and (3{t) = tt. Therefore the calculus of the different elements defined in section 2 give 



, , / / X , t sin2t 

a(t) = - Pisjds = - 2 + ' 



1-e"^ 



a{T) = ao = Y^^^ j /3(s)e-'^(^)ds « 6.752 
- (137^7 X /^(^)e-"^^^^^ « 4.0912 



In fig.l we have drawn the graph of the function t 7(i)i 7(0 = o-a + f3{s)e~°'^'^^ ds 

L ./o 

8.5 



,a{t) 




Figl. The periodic function 7(i) 
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The fig.2 indicates the graphs t yoo(t) and t^yoo{t + T) where yoo(i) = 1{T) + b{s)e °^*^ds 

L Ja 

60 



,a(t) 



40 



30 



10 



-10 




2 4 



8 10 12 14 16 18 20 

t 



Fig2. The asymptotic function yoo(i) 



The fig.2 fairly shows the drift property of yoo{t) just as the function b{t). 
The asymptotic behavior represented in fig. 3 shows us the very fast convergence of the approximation 
Zn{t) = Vooit) + n"f{t) to the exact solution at the point t + nT i.e. yn{t) = y{t + nT), y{t) being the solution 

of the Cauchy problem (6.1) and given by y{t) = \y{0) + / b{s)e-''''''Us^ f^"'-*^ 



So for n = 1 the two graphs coincide exactly for t > 5! 
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Fig3. Asymptotic behavior, n=l 

Finally in fig. 4 we have put the graphs of the functions yn{t) and + tt) with n = 5 and here we also note 
the drift property for the solution of (2.1) taking initial value yo = 1 at t = 0. 
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Fig4. Drift property in asymptotic behavior 
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